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Abstract 



A new C-integrable limit of the second harmonic generation equations is found. 
The corresponding general solution is given in an explicit form. Connection of this 
problem with the modified Liouville equation is discussed. 

d, 

Second harmonic generation (SHG) is a well-known method of experimental investiga- 
tions in nonlinear optics. Under idealized condition of only one space dimension this process 
is described by the equations (see, e.g., 

Qix = ~ 2 <l2ql, q2r = q\, (1) 

P 

>• ■ where differentiation with respect to characteristic variables x an d r is indicated by sub 
scripts, the asterisk denotes complex conjugation, and q\ and q 2 are slowly varying com- 
plex electric field amplitudes of two electromagnetic waves with carrier frequencies u>i and 
0J2 = 2oq, respectively. The SHG equations (JJ) are shown to be "S'-integrable" |2| what per- 
mitted one to find some particular solutions describing transformation of waves. However, 
application of the inverse scattering transform method to the Cauchy problem meets some 
peculiar difficulties and is not fully developed yet. Therefore, other approaches to solving 
the SHG equations are of considerable interest. 

As it was noticed long ago ||, in the case of pure amplitude modulation, when q\ and q% 
are real variables, the SHG equations can be reduces to the C-integrable Liouville equation, 
so that the solution can be expressed explicitly in terms of two arbitrary functions. Some 
particular examples of such solutions were studied in J3J, and the Liouville solution was 
applied to the so-called restricted Cauchy problem typical for experimental situation in []5J. 

Here we want to note that the SHG equations (HD can be approximated under certain 
condition by the system which is C-integrable and its solution is given in an explicit form in 
terms of four arbitrary functions. This limit is in a sense opposite to the discussed previously 
pure modulation case. (Note, however, that a pure modulation case is an exact reduction 
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of the SHG equations, whereas here we consider some approximation to these equations.) 
The proposed approximation is based on the following observation. Elimination of q 2 from 
Eqs. (0) leads to the equation 

qi, X rQi ~ 9i, x 9i,t = -2 (<?i<?i*) 2 , (2) 

where the left hand side is quadratic and the right hand side is of fourth degree in the 
amplitude q\. Hence, one may suggest that in the limit of small enough \qx\ we can neglect 
the right hand side, if derivatives of q\ in the left hand side are not too small. To formulate 
the corresponding criterion, it is convenient to pass to real variables introduced in ||. We 
represent q\ in the form 

qi = (VQ/2) exp(z0/2), Q > 0, (3) 

so that q 2 is given by 

q 2 = ~\ (\nQ) x + i<P x exp(i0), (4) 
and Eq. (||) reduces to the system 

(ing) XT -0 x T = -g, (Q0 T ) x = o. (5) 

Introducing the variables 

s = l/Q, u = 4> T , v = (\ns) x , w = (f) x (6) 
we rewrite the system (||) in the form || 

(s) x — sv, u x = w T = vu, v T = l/s — uw. (7) 

If we put here u = 0, we obtain at once the Liouville equation (In s) XT = 1/s which application 
to SHG was considered in 0, [5|. On the other hand, if 1/s <g \uw\, or 

Q<\4> T 4> X \, (8) 

we arrive at the system 

v T = —uw, w T = uv, u x = uv, (9) 
whereas s can be found from the equation (ln(s/w)) x = which follows from Eqs. (0), i.e. 

s = K(t)u, (10) 

where K{r) is an arbitrary function. 

To find the general solution of the system (Q), we eliminate v and w, 

v = — , w = — (^) =--^ + ^1, 11 
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so that w T can be presented in the form 



UK 



1 [U TT - Uu 2 T /u) 



a 



According to Eq. (PI), this must be equal to u x , hence we obtain 



3f„,2 



u 2 T /u) 



or 



u 



3^ 

2 u 



-l(« 2 )- 



1„.3 



(12) 



where /(r) is an arbitrary function. Generally speaking, this equation cannot be solved for 
a given /(r) in an explicit form. However, in our case /(r) is an arbitrary function and it 
can be replaced by another arbitrary function F(t) related with f(r) in the following way. 
We make a substitution 



u = F>(t)U(0, Z = F(t), 



(13) 



so that Eq. (O) reduces to 



U, 



\ul/u) + \u 



1t T 3 
2 



2V t'"' ' 2 ~ (F'(r)) 2 



/(r) 



3 (F"(t) 
2 V /'(r) 



F"'(t) 
F'(r) 



U. 



Thus, if f{r) is expressed in terms of F(t) by the equation 



f(r) 



F"'(t) 3 /F"fr 



F'(r) 



fir) 



(14) 



where the curly bracket denotes the Schwarzian derivative |7j, then U(£) is determined by 
the equation 



0. 



(15) 



This equation is solved by elementary methods to give 

U(0= [l/G + \(Z-H) 2 ]-\ 

where G and H are integration constants (arbitrary functions of x)- Making use of Eqs. ([13 
and (|ll|) , we arrive at the general solution of the system @ : 



u(x,r) 
w(x,r) 



AF'G 



A + G 2 (F-H) 2 ' 

—GG'(F - H) 2 + 2G 2 E'{F - H) + AG'/G 
A + G 2 (F-H) 2 ' 
\G Z H'(F - H) 2 + AG\F - H) - 2GH' 
4 + G 2 (F-H) 2 ' 



(16) 
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where F(r), G(x), H(x) are arbitrary functions. Together with Eq. (0), these formulas 
give full description of the SHG problem in the limit (§). 

^From mathematical point of view, the fact of C-integrability of the system (^) is not 
trivial and finds its explanation in possibility to reduce this system to the so-called modified 
Liouville equation (see Appendix). 

One may hope that this solution will permit one to give analytic description of SHG 
processes for concrete examples in a way similar to that used in ||. 

We thank A.V. Zhiber for indication of Ref. || and H. Steudel for useful remarks on 
the manuscript. AMK is grateful to DFG (grant 436 RUS 113/89/2 (R,S)) and INTAS 
(grant 96-0339) for partial support. MVP is grateful to RFBR (grants 00-01-00210 and 
00-01-00366). 

Appendix 

Let us show that the system (j^) can be reduced to the modified Liouville equation first 
introduced by E. Vessiot || (see also the paper [[|). 

We notice that the last two equations @ yield a first integral 

v 2 + w 2 = W( X ))\ (A.l) 
where ip(x) > is an arbitrary function. Then transition to polar coordinates 

v = i>'(x) cos 9, w = ip'(x) srn 0j (A.2) 
transforms the last two equations @ to 

u = 9 T , (A.3) 

and the first equation (^) takes the form 9 TX = ip (x) cos 9 ■ 9 T or 

9 rC = cos9-9 T , (A.4) 

where ( = ip(x)- (This equation can be integrated once and reduced to a linear second order 
equation; see Eq. (|A.20|) below.) Introducing a new variable 

(ft = \nu = \n9 T , (A.5) 

we obtain 

= ®t<;/9 t = cos 9, 0£ T = — sin 9 ■ 9 T , 

and finally 

rC = -e^l-tPl (A.6) 

which is just the modified Liouville equation [|). (Amplitude-modulated reduction of de- 
generate two-photon propagation equations discussed in [|10j can be also reduced to this 
equation.) 
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This observation permits us to write down the solution of Eq. ( |A.6[) . Indeed, the general 
solution ([16D gives 

(^(x)) 2 = v* + w 2 = (G> /Gf + \{GH')\ 

Hence, if we choose ip'(x) = 1j so that ( = x an d H{x) * s connected with G(x) by the 
equation 



m 2 = -^[g* - 



then the solution of Eq. 



is given by 
6 = hiu = In 



AF'G 



4 + G 2 (F — H) 2 ' 



(A.7) 



(A. 



Note also that another system (arisen first in a physical problem about amplification of 
two counter propagating light beams by transverse flow of an amplifying medium) 



a T = ac, b T = —be, c y 



(a + b)c, 



whose general solution was found long ago [p] , can be reduced by substitutions 

u = c, w = b — a, v = —(a + b), 

to the system 

v T = uw, w T = uv, u x = uv, (A. 9) 

which differs from Eqs. (^) only by the sign in the right hand side of the first equation, and 
is connected with another forms of the modified Liouville equation 



depending on the choice of the sign in the first integral 

v 2 -w 2 = T^\x)f- 
General solution of Eqs. ( |A.9 ) is given by the formulas 

AF'G 



(A.10) 



(A.H) 



w(x,r) 



4 — G 2 (F — H) 2 ' 

GG'{F - H) 2 - 2G 2 H'(F - H) + AG'/G 
A-G 2 {F-H) 2 ' 

-\G Z H'(F - H) 2 + AG'(F - H) - 2GH' 
4-G 2 {F-H) 2 ' 



(A.12) 



Then we have 



v 2 -w 2 = G' 2 /G 2 -\G 2 H'\ 
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and 



is the solution Eqs. ( |A.10| ) provided H satisfies the equation 

H' 2 = A\G' 2 ± G 2 ]/G 4 (A. 14) 

with a proper choice of the sign. 

As was noted in ||, solutions of the modified Liouville equations ( |A.6| ) are connected 
with solutions of the Liouville equation 

z T< = e z (A.15) 

by differential substitutions 



3 = + ln(0 f + ^/$±l), ( A - 16 ) 
which can be written in terms of hyperbolic functions as follows, 

z = (p + Arsinh0£ or z = <fi + Arcosh0^. (A. 17) 

These substitutions map the general solutions ( A.13|jA~T4"D of the modified Liouville equations 
( |A.10p into the general solution 

WW+QKif' (A ' 18) 

of the Liouville equation (|A.15|) , where 

Q(0 = -2/G(0-H((). (A.19) 



Substitutions (|A.17|) can be reduced to a second order linear equation for the variable y 
defined by y^/y = \ exp(z — (p) which is considered as a function only of £: 

The above consideration gives one solution of this second order equation in the form 

y = exp n J exp(z - (j>)d{J . (A.21) 

At last we notice that Eq. ( |A.6|) does not admit real substitution ( |A.17| ) (but only complex 
one) which transforms this equation into Liouville equation. Nevertheless, if we write the 
solution of the modified Liouville equation 



rC = -e*y/ 1 - e^l (A.22) 

in the form 

2F'(r)G(() 

0_ln (F(r)-i/(C)) 2 + ^(C)' (A ' 23) 

where 

H' 2 =G 2 -e 2 G' 2 , 

then in the limit e: — ^ it goes directly to the solution of the Liouville equation T( * = — exp 0. 
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